Ultrafast spherulitic crystal growth as a stress- induced phenomenon 
specific of fragile glass-formers 
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We propose a model for the abrupt emergence, below temperatures close to the glass transition, of the 
ultra-fast (GC) steady mode of spherulitic crystal growth in deeply undcrcoolcd liquids. We interpret this 
phenomenon as controlled by the interplay between the generation of stresses by crystallization and their 
partial release by flow in the surrounding amorphous visco-elastic matrix. Our model is consistent with both 
the observed ratios (~ 10 4 ) of fast-to-slow velocities and the fact that fast growth emerges close to the glass 
transition. It leads us to conclude that the existence of a fast growth regime requires both (i) a high fragility 
of the glassformer; (ii) the fine sub-structure specific of spherulites. It finally predicts that the transition is 
hysteretic, thus allowing for an independent experimental test. 



Upon approaching T g , the rate of crystal growth from 
deeply supercooled glassformcrs decreases to very small 
values - in the 10 -12 m/s range - as appears consistent 
with the dramatic slowing down of molecular motions. 
Growth is hence expected to be nearly arrested in the 
glass phase. This expectation, however, is challenged by 
a surprising phenomenon identified by Oguni and cowork- 
ers in several fragile materials 1 ' 2 and further investigated, 
more recently, by the groups of Yu and Ediger at Madi- 
son 3 . Namely below some temperature T t slightly larger 
than T g , an ultrafast crystal growth mode (often coined 
GC, for Glass-Crystal) suddendly emerges. Near T t the 
ratio of the fast over normal (slow) front velocities, Uf 
and u s , is huge, typically of order 10 4 . 

Up to now this phenomenon has been observed in 11 
one-component molecular glass- formers 3 , which all have 
high fragility indices, that is, exhibit the so-called decou- 
pling phenomenon 4 ' 5 . Namely, below ~ 1.2T g the viscos- 
ity r\ and translational diffusion coefficient D of the SC 
liquid no longer obey the Stokes- Einstein relation but in- 
stead, D ~ r]~*, with < £ < 1. Ediger et al 6 have 
shown that, as T varies, the growth velocity of the slow 
(normal) mode scales, not as ?y _1 , but as D, a result 
consistent with the idea that crystallization kinetics is 
controlled by local processes 7 . 

Moreover, materials exhibiting fast growth share the 
following common features: molecules are strongly 
anisotropic (nearly planar); single-crystals grown at 
small undercoolings arc highly facetted - which indicates 
slow interfacial kinetics; at large undercoolings, crystal 
growth is spherulitic 8 ' 9 with, in the usual observation 
domain (spherulite radii ~ 10's to 100's of microns), a 
quasi-constant front velocity. 

To this day, the paradoxical emergence of this fast 
growth mode near T g largely remains a mystery A no- 
ticeable proposition formulated by Tanaka 10 is that be- 
low T g , volume contraction upon crystallization creates 
in the surrounding glassy matrix a negative pressure, i.e. 
free- volume; this would increase particle mobility close 



to the interface, hence accelerate crystallization. Konishi 
and Tanaka 11 have shown that the values of the veloc- 
ity "jump" at T t between fast and slow growth for two 
polymorphs of salol indeed correlate with their density 
changes upon crystallization. The Madison group, how- 
ever, have later objected that such is not the case for the 
polymorphs of another material (ROY) 12 . They more- 
over point out that a net accumulation of free- volume at 
the growing interface would result in accelerated growth, 
contrary to observations. Whether GC growth is a stress- 
induced phenomenon thus remains an open issue. 

We will here uphold that such is indeed the case, pro- 
vided that Tanaka's ideas are formulated in a slightly 
more general way. Indeed, it seems physically undis- 
putable that: (i) volume contraction upon crystalliza- 
tion creates stresses in the surrounding matrix and at 
the interface; (ii) these stresses modify interfacial kinet- 
ics. Now, we will argue that: (i) contraction does not re- 
sult in net free- volume creation at the interface; (ii) how- 
ever, in the absence of an average negative pressure, due 
to structural disorder, deviatoric interfacial stresses and 
pressure fluctuations are sufficient to accelerate growth. 

Besides, we note that GC growth appears near T g , in 
a region where, as the relaxation time r Q increases, the 
mechanical response of the amorphous matrix rapidly 
crosses over from predominantly viscous (above T g ) 
to predominantly elastic. Its mechanical response is 
hence visco-elastic. The matrix rheology determines to 
which extent the elastic Eshelby stresses generated by 
spherulitic growth are relieved by the dissipative flow. 
The interfacial stress is thus controlled by the competi- 
tion between the generation of stresses by growth and 
their relaxation by flow in the matrix. 

In the following, we present a schematic model, which 
describes this interplay in a simple way, using the lin- 
ear Maxwell rheology to model the matrix mechanical 
response. The model is developed in two stages. In its 
first "idealized" version, we represent the spherulite as 
an homogeneous, elastically isotropic, sphere. We then 
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predict at low enough temperatures the coexistence of 
two growth modes: 

- a slow mode, which continously extrapolates from high 
temperatures, and in which interfacial stresses are almost 
fully relaxed by viscous flow in the matrix; 

- a fast one, in which, on the contrary, viscous stress re- 
laxation is inefficient, so that stresses are essentially those 
produced by an Eshelby inclusion in an elastic matrix. 
This model satisfactorily accounts for the order of mag- 
nitude (typically ~ 10 4 ) of the ratio between the corre- 
sponding velocities, Uf and u s , yet also predicts that a 
"fast spherulite" would commute to slow growth past a 
T-dependent maximum radius which turns out to take 
unphysically small values (in the 100 nm range) around 
both T g and T t , in contradiction with observations. 

This leads us to question the above simplistic descrip- 
tion of a spherulite. Indeed, as documented in detail 
by Faivre and coworkers 9 , spherulites are compact ob- 
jects which appear spherical on the optical scale, but 
are constituted, at smaller scales, of a space-filling accu- 
mulation of "sheaves" (of micrometric diameters), them- 
selves constituted of anisotropic single crystals (typical 
size ~ 100 nm), growing intermittently via small- and 
large-angle branching. Therefore, although spherulites of 
more than micrometric sizes can be viewed on a coarse- 
grained scale as spherical isotropic elastic inclusions in an 
amorphous matrix, crystallization actually proceeds via 
the intermittent growth of fiber-like units (the sheaves). 
This leads us to propose a second version of the model, 
which takes into account the fact that the growing units 
have a constant head radius comparable with a fiber di- 
ameter. This "rugged spherulite" model now predicts 
that fast growth emerges at low temperature via a hys- 
teretic, saddle-node bifurcation, the existence of which is 
a direct consequence of the decoupling phenomenon. It 
accounts not only for the ratio of growth velocities, but 
also for the observability of fast growth up to the vicin- 
ity of T g . Our interpretation of the GC-vs- normal growth 
phenomenon in terms of elastic vs viscous stress control 
moreover explains an important observation made by the 
Madison group, namely that the fast growth velocity is 
comparable with that of solid-solid transformations be- 
tween polymorphs. 

In Section I we derive the equations which relate the in- 
terfacial stress and the front velocity for an elastic sphere 
growing from a visco-elastic matrix. In Section II we 
show that such an idealized spherulite may present fast 
growth, but this only transiently, and up to maximal radii 
which are much too small to account for experimental 
observations. In Section III we describe our modified, 
"rugged spherulite" , model and show that it accounts for 
the emergence of steady fast growth at low temperatures 
via the crossing of a sub-critical, hysteretic, bifurcation. 
In the discussion (Section IV) , after summarizing our ar- 
gument, we conclude that two physical ingredients are 
indispensable for the emergence GC growth, namely a 
high fragility of the supercooled liquid and a fine sub- 
structure of the growing crystal phase. 



I. GROWTH OF A SOLID SPHERE FROM A 
VISCO-ELASTIC MEDIUM 

In the temperature range where GC growth occurs, 
close to the glass transition, crystal front velocities are 
at most of the order of 10~ 7 m/s. The temperature rise 
ST at the front due to latent heat production is of order 
VCR/ K m 13 , with V the front velocity, R the spherulite 
radius, C the latent heat per unit volume, and K m the 
thermal conductivity of the amorphous matrix. With 
C ~ 10 8 ,I/m 3 , K m ~ 10- 1 J/(m.s.K), and for typical R 
values < 200/^m 2 ' 12 , we get ST < 10~ 2 K. We can there- 
fore consider growth to occur in a completely isothermal 
matrix. 

In this section, we assume that the crystal is an 
isotropic elastic solid with a perfect spherical shape, 
growing from a mechanically isotropic, linear visco- 
elastic, amorphous matrix. Under the ambient pressure 
POj the matrix and crystal have equilibrium densities po 
and p C Q = po + Spo, with 5po > 0. That is we assume 
that the strain associated with crystallization is a pure 
volume contraction. 



A. Interfacial stress dynamics 

From a mechanical standpoint, we are dealing with an 
Eshclby-like problem 14 , as is clear if we consider, for a 
moment, the case when the matrix is a perfectly elastic 
solid. The stress field r can be computed as follows: at 
any given instant t, a spherical volume Vq of the matrix 
has transformed into a (crystalline) state characterized 
by a strain Iq = — 1- Mechanical equilibrium then 
demands that V -r = inside each phase (i = c, m) with: 

f l = 2 pi ti + A, trei T (1) 

and 

% = \ (vNZ, + VuJ ) (2) 

where is the elastic strain relative to the reference state 
of each phase (equilibrium under pressure p), and Hi the 
corresponding displacement. These equations must be 
supplemented with two interfacial conditions: (i) me- 
chanical equilibrium, i.e. continuity of normal stresses; 
(ii) mass conservation on the moving front. 

Since we are in spherical geometry under pure contrac- 
tion, it is clear that: 

Ui = u^r)^ (3) 
Mechanical equilibrium then entails that: 

in the matrix (r > R) : u m (r) = A m /r 2 (4) 
in the crystal (r < R) : u c (r) = A c r (5) 

It immediately follows that V • u m = everywhere in 
the matrix, in particular at the interface. This leads to 



3 



a classical 15 but counterintuitive result: volume contrac- 
tion upon crystallization does not give rise to any excess 
pressure in the surrounding matrix; it creates a purely 
deviatoric stress field. 

The above discussion immediately carries over to the 
case of our concern of a linearly visco-elastic matrix, as 
the sole difference is that the Lame coefficients A m and 
p m now become frequency dependent. The space de- 
pendence of the displacement fields is still given by (4) 
and (5) (with time-dependent coefficients A m and A c ). 
Hence, even in the presence of a dissipative mechani- 
cal response, the stresses remain purely deviatoric in the 
amorphous matrix. It follows that the viscous flow, if 
any, is also divergence free. Finally, no free-volume is 
created. 

In the crystal (r < R), the stress field (a pure pressure) 
is given by: 



r c (r,t) = (3K c A c - Po )l 



(6) 



with K c = A c +|/x c the crystal bulk modulus. The crystal 
is under strain A c 1, i.e. grows at density p c o + 5p c such 



that 



-4, 



ISpc 

3 Pco 



(7) 



In the matrix (r > R), the strain field (purely devia- 
toric) is: 

%n(r, t) = —^-(-2e r e r + e g e e + e<f,e$) (8) 

Since we assume a linear Maxwell rheology the stress in 
the matrix is necessarily of the form: 



7 ~ l(T 2 N/m, K c > 10 10 Pa, and R > lpm, this is at 
most of order 10~ 6 , hence completely negligible. 
(ii) mass conservation across the interface moving at ve- 
locity R: 



Po {u m (R)-R) = Pc (uc(R)-R) 



(13) 



with, for any r, u m ( r ) = 4^" anc ^ 
the above equation also writes: 



A r r, so that 



A m 



PO-PcO-SPc\R , PcQ+bPc 



Po 



R 



A c 



Po 



(14) 



The relative density change upon crystallization e = 
Spo/po is at most a few percent. Since we know from 
equation (7) that A c is 0(e) we neglect the quadratic 
term -^-A C {R) in the above equation, whence: 



Am / 



R 

■e + 3A C ) -+A C 
H 



(15) 



Together with equations (11) and (12) this now provides 
an equation relating the interfacial stress to both the 
growth velocity and the spherulite radius: 



(16) 



where we have defined an equivalent elastic modulus p 
via: 



1 1 



2 p 2p m 3K C 



(17) 



T m (f, t) = r(r, t)(-2e r e r + e e e e + e^e^) -p 1 (9) 
where t verifies: 



0-kr{r,t) 



n r^l^ (10) 



2p m dt 2rj 



- 3 dt 



with rj the viscosity of the amorphous phase. This solves 
at any point r as r(r, t) = O -1 * A m wherefrom the 
interfacial radial (tensile) stress S(i) = — 2r{R{t),t) ver- 
ifies O* (S R 3 ) = -2A m , i.e.: 



2p m dt 



t dA n 
' dt 



(11) 



Mechanical equilibrium in the bulks of both phases 
must now be complemented by two interfacial conditions, 
namely, for our growth problem: 
(i) continuity of normal stresses: 



B. Interfacial kinetics 

Equation (16), which expresses both mechanical equi- 
librium and mass conservation conditions, must now be 
supplemented by the specification of interfacial kinetics. 
For a simple liquid, the crystallization process can be 
described as a series of activated jumps over barriers, 
of typical height E, separating the amorphous and crys- 
talline molecular configurations. At very large under- 
cooling, the probability of back jumps from the crystal 
into the matrix is negligible so that one may write the 
front velocity (in the absence of any interfacial stress) as 

E 

~ ljq a e k s T , with uiq a typical molecular frequency and 
a an atomic length. In the fragile glassformers of interest 
here, molecular transport involves a wide, T-dependcnt, 
distribution p(E; T) of barrier heights, so that one should 
schematically write: 



3K C A C = £ 



(12) 



Note that we have omitted here the contribution of the 
interfacial tension 7. Indeed, it introduces a correction 
to the internal strain of the crystal of order ■ Taking 



R = V (T) ~u a / dEp(E; T) e "b 



(18) 



Vq(T) is the growth velocity that should be measured in 
the absence of stress effects. 
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FIG. 1. Phase portrait of the dynamical system (22) in the (r,a) plane, for b = 3.5, 5, and 10 (from left to right). 



As discussed above, although the growth of an ideal- 
ized spherulite produces no free- volume, it does generate 
finite shear stresses in the interfacial region. Like glasses, 
deeply undercoolcd liquids close above T g exhibit highly 
multistable free-energy landscapes. As elaborated long 
ago by Eyring, stresses bias such landscapes: a series of 
works 16-19 has indeed now clearly shown that this effect 
controls the plastic deformation and rhcology of amor- 
phous solids. When an amorphous system is subjected 
to a given external shear stress field, the barrier limit- 
ing a particular local reconfiguration may increase or de- 
crease depending on the relative orientation of molecular 
displacements with respect to the principal axes of the 
stress. 

In the presence of interfacial shear stresses S, a barrier 
of height E is thus shifted into E + SE(T,) and follow- 
ing Eyring, we assume that 5E = aT,. Given the high 
structural disorder of the amorphous phase, we can rea- 
sonably consider that the a's have random signs. For the 
sake of simplicity, we furthermore assume that they are 
sharply peaked around ±ao, with ao a typical "activation 
volume" . The growth front velocity now becomes: 



R = V (T) cosh 



(19) 



where Vq(T) is given by equation (18). Our growth prob- 
lem is now fully specified by equations (16) and (19). 



II. GROWTH OF AN IDEALIZED SPHERULITE 

In this section, we analyze the solutions of the above 
growth model. That is, we solve equations (16) and (19) 
starting from initial conditions corresponding to an al- 
ready formed idealized spherulite well beyond the nucle- 
ation regime. Since we are primarily interested in the 
transition between normal and GC growth, we restrict 
our attention to temperatures close to T g . 

We first rewrite the growth equations in a dimension- 
less form, by introducing the following units for, respec- 



tively stress, length, and time: 



IT 



x= v -v 



and 



/' 



/' 



(20) 



Note that, in contrast with the stress scale Eo, both A 
and r depend on temperature via the rapidly varying 
quantities r/ and Vo- Denoting 



and 



R 



equations (16) and (19) now become: 

3r ,„ , 
a + a = — (1 - a) 
r 

f = cosh(6 a) 



(21) 



(22a) 
(22b) 



The reduced dynamics is thus controlled by a single pa- 
rameter 



b = 



(23) 



which measures the sensitivity of interfacial kinetics to 
tensile stress. 

The behaviour of the solutions of the dynamical sys- 
tem (22) is best summed up by plotting its phase portrait 
in the (r, a) plane, as shown on Fig. 1 for three values of 
b. The (j = condition, which is verified when: 



r*(a) = 3- — -cosh(6o-) 



(24) 



defines a curve, denoted C*, in the (r, a) plane (dashed 
lines on Fig 1). On the left (r < r*(a)) of C* , the inter- 
facial stress increases as growth proceeds. On its right 
(r > r*(cr)), a decreases until it approaches asymptoti- 
cally the low-cr, high-r part of C* . 

Starting from a small initial r value, the stress a first 
increases very sharply. For small 6's, it quickly drops 
back to very small values after reaching its maximum (at 
the crossing point with C*). 

As b increases beyond b c ~ 4, C* develops an increas- 
ingly marked reentrance with a flat upper branch a ~ 1 
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extending up to a value r max that grows exponentially 
with b. The trajectories stay close to this high-er branch 
up to radii > r max , then drop sharply to reach the flat 
(7 <C 1 branch. Note that for b 3> b c (see Fig lc for 
b = 10), as most trajectories originating from the left of 
C* (see curves 1 to 5 on the figure) come closely together 
onto the common high-cr plateau, they also drop at about 
the same value r sw i tc h- 

Since [see Eq. (22b)] a measures the logarithm of the 
growth velocity, the high-cr plateau corresponds to a tran- 
sient fast growth regime extending over a range of radii 
which, although finite, increases exponentially with b. It 
now appears natural to check whether we can interpret 
the upper and lower branches of C* as corresponding to 
the GC and normal spherulitic growth regimes respec- 
tively. For this purpose, we must be able to account 
first of all for the two main experimental features: (i) 
the order of magnitude of the ratio of the fast to normal 
growth velocities in the vicinity of T g ; (ii) the existence 
of an upper temperature T t ^ T g above which compact 
GC growth is not observed. 

Comparisons with experiments are here performed 
on the especially well documented case of OTP (o- 
tcrphcnyl). For this material, at T g = 243K 20 , where 
the viscosity rj = 10 11 Pa.s, the contraction parameter 
e « 4.8% 21 ; using the longitudinal wave speed cj = 
2550m/s 22 and the ratio q/q ~ 2 23 we estimate a shear 
modulus value /i ~ 1.7GPa; the molar mass is M = 
230g/mol; and the specific volume v = 0.89cm 3 /g 24 . 

(i) Velocity ratio. In our idealized spherulite model, 
the normal growth velocity (cr -C 1) is u s ~ Vq. Fast 
growth (cr ~ 1) occurs under an interfacial tension ~ 
Eo = -4f±, corresponding to a shear stress ~ So/2 = 

~ 3.1CP 2 /i - that is, a sizeable fraction of the yield 
stress. The growth velocity Uf ~ Vbcoshfr. From the 
measured value Uf/u s ~ 10 4 , we deduce that b 10. 

A first assessment of whether fast growth may result 
from the stress-induced acceleration of interfacial kinetics 
consists in checking that these parameters correspond to 
a reasonable value of the activation volume a^. With 
the help of Eq. (23), we find «o ~ 0.3 nm 3 , precisely of 
the order of the molecular volume v m ~ 0.34 nm 3 , which 
indeed brings valuable support to this idea. 

(ii) Switching from GC to normal growth. Our ide- 
alized spherulite model predicts that GC growth (when 
initial conditions allow for its occurrence) cannot be a 
permanent regime, but must switch to normal growth 
when the reduced radius reaches an upper limit r SW i tc h- 
Can we reconcile this prediction with the experimental 
observation of a temperature T t above which compact 
GC growth ceases to be observed? Let us note that 
the published data arc concerned with spherulite radii 
of a few hundreds of micrometers at most, hence do not 
necessarily contradict the possibility that GC growth be 
transient. 

The upper physical radius at which GC growth is pre- 



dicted to be observable is 

-Rswitch = A r switc h(&) (25) 

where [see Eq. (20)] A = r)Vo/(i. As shown by the Madi- 
son group 25 materials where GC growth is observed have 
high fragilities and, as such, exhibit a decoupling between 
translational diffusion and viscosity 5,26 ' 27 . Namely, be- 
low about 1.2 T g the translational diffusion coefficient D 
no longer scales as ry -1 , but as f]~Z with < £ < 1. As 
emphasized in 25 , these materials have comparable values 
of £ of order 0.7. Moreover, Ediger et al. 6 have shown 
that the normal growth velocity Vq scales, not as rf~ , 
but as D - a result consistent with the idea that crystal- 
lization kinetics is controlled by local processes. In the 
narrow vicinity of T g we are interested in, the equivalent 
shear modulus fi is quasi-constant as compared with the 
fast- varying quantity n Vq and our lengthscale A scales 
with T as: 

X^^D^r] 1 -^ (26) 

As temperature increases, b [see Eq. (23)], hence 
''switch (b), decreases and so does A. At a fixed tempera- 
ture GC growth would thus be observable up to a max- 
imum radius i? S witch, which decreases with T. We must 
now examine whether our idealized model is compatible 
with the observation of sphcrulitcs of radii ~ a few 100/im 
at T g . For b = 10, we find r sw i tc h — 3.10 3 ; from 20 we es- 
timate for OTP Vq ~ 10~ 12 m/s; with the above values 
of rj and /i, this yields: i? S witch ~ 0.2^m. Clearly, the 
idealized spherulite model, as such, fails to account for 
experimental observations. 

Does this failure invalidate our primary assumption 
that GC is a stress-induced phenomenon? Or could it 
be attributed to shortcomings in our extremely idealized 
representation of a spherulite? 

It is useful at this stage to return to Eq. (16) and note 
that when stress relaxation in the amorphous matrix can 
be neglected, that is for fi/rj = 0, the interfacial stress 
asymptotes to the stationary value £o = 4e/j,/3 28 , which 
is precisely the stress level at which transient GC growth 
occurs in the above description. This means that, in the 
stress-controlled model, GC growth is too fast for stress 
relaxation in the supercooled matrix to bc efficient and 
occurs at a velocity close to that of a crystal growing 
from a solid matrix. This qualitative result must be put 
in regard with an important remark made by the Madi- 
son group 12 : namely, they point out that the order of 
magnitude of GC growth velocities compares with that 
of solid-solid transformations between polymorphs of the 
same material. 

In our opinion, this, together with the internal con- 
sistency of our prediction for the velocity ratio, makes a 
strong case for retaining the stress-controlled assumption 
and reconsidering the details of our model. 
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III. GROWTH OF A "RUGGED SPHERULITE" 

The representation we have used so far of a sphcrulitc 
as a perfectly spherical homogeneous solid was drasti- 
cally over-simplified. Indeed, these objects result from 
the growth of anisotropic micro-crystals which, as they 
multiply via intermittent small- and large-angle branch- 
ing, form a space-filling quasi-spherical polycrystallinc 
structure. So, as seen on Fig. 6b of reference 9 , the crystal- 
matrix interface is rugged on a scale in the ^100 nm 
range, fixed by the micro-crystal growth habit, and thus 
independent of the average sphcrulite radius R. Clearly, 
ahead of such an interface, the stress field must be mod- 
ulated on the same scale. 

We must therefore reformulate our growth model to 
take into account the presence of this interfacial rugged- 
ness. Our model boils down to two equations: (i) the 
kinetic equation (19) which relates the growth velocity 
to the interfacial stress; (ii) the mechanical equation (16) 
which describes the evolution of interfacial stress under 
the combined effects of growth and viscous flow in the 
matrix. The kinetic equation, which is local, hence in- 
dependent of the inclusion shape, clearly remains un- 
changed. The mechanical equation, however, which was 
derived in Section I for the special case of a perfect 
sphere, must now be reconsidered. 

For this purpose, we first reinterpret equation (16) in 
the simple, fully elastic, limit (fi/rj = 0). We denote C(t) 
the domain occupied by the crystal at time t, and I(t) 
the corresponding interfaces. Let us assume that at time 
t\, the stress field in the system satisfies the conditions of 
mechanical equilibrium (V-r = everywhere; continuity 
of normal stress on I(ti)). The change of the stress field 
due to spherulitic growth over the interval [£1,^2] is due 
to the amorphous-to-crystal (A — > C) transformation of 
the domain extending between C(ti) and C(^)j thanks 
to linear superposition, this is equivalent to the following 
sequence: 

1) C —> A transformation of C(ti) 

2) A —> C transformation of Cfa) 

In other words, inside the matrix, at any point: 

r(f,i 2 ) -?(r,ti) = r Es h.(r, t 2 ) ~T Es h.{r,ti) (27) 

where TE S h.{r 7 t) is the Eshclby field 14 ' 15 due to the A — ► 
C transformation of an inclusion C(t). 

Consider now a point r(t), which lies on the interface 
at all times, and is constrained to move along the local 
normal vector (df/dt = (V ■ n)n). Using (27), the inter- 
facial stress at this running point, r(f(t),t), is found to 
verify: 

df(r(t),t) (f . = v df Esh .(r(t),t) 
—— = {V . n)n .V( T -r Esh .) + 

(28) 

Coming back to the case of a growing sphere, for which 
7 = Esh.( , r(t) 5 1) = Cst, let us rewrite equation (16), which 



governs the dynamics of interfacial stress, as: 

E + ^£ = ^(£ -I]) (29) 
77 R 

The second term on the l.h.s of this equation results from 
our assumed Maxwell rheology in the matrix. Compar- 
ison with Eq. (28) then shows that the factor R/R in 
the r.h.s. must be understood as the ratio of the normal 
growth velocity to the characteric length scale of stress 
gradients ahead of the growing front. 

Clearly, given the complexity of spherulitic growth , a 
detailed description of the rugged front geometry as well 
as of the map of local normal velocities remains for the 
moment out of reach. Let us observe, however, that, as 
mentioned above, a spherulite can be viewed as a com- 
pact sphere of radius R constituted of space-filling sub- 
structures of typical lateral dimension D ~ lOOnm (we 
refer here to the scale denoted e in Rcf. 9 ). We evalu- 
ate in the appendix the stress field near the surface of a 
spherical inclusion decorated with bulges of typical size 
D « R. We find that: (i) the tensile stress at the 
interface remains nearly equal to its value So for the 
ideal sphere; (ii) the dominant length scale controlling 
the stress gradient is the typical bulge size D. 

On the basis of these remarks, we hypothesize that 
the typical interfacial stress is governed by an equation 
similar to (29), yet in which the factor R/R is replaced 
by V„/D = R/D, where D is a constant characteristic of 
the spherulitic micro-structure: 

S + ^S = ^(£o-£) ( 3 °) 
77 D 

while the kinetic equation (19) remains unchanged. In 
the reduced units defined in (20), and after elimination 
of R, the dynamical problem reduces to: 

3A 

a + a = — (1 - a) cosh(kr) (31) 

which is similar to Eq.(22a), with the difference that the 
dynamical variable r is now replaced by the constant 

d = D/X (32) 

In contrast with the ideal spherical growth model 
[Eq. (22)], the dynamical problem (31) presents fixed 
points, which are defined by the condition a = 0, or 
cquivalently by: 

d = 3- — — cosh(6cr) (33) 
a 

The loci of these fixed points in the (d, a) plane, a* (d) , 
are the same as the dashed curves shown on Fig. 1 and 
are reproduced on Fig. 2. The fixed points are stable 
when da* /dd < 0. When b > b c ~ 4, as shown in 
Section II, the curves a* develop a reentrant branch, 
which is thus unstable. In this case, the nearly flat high- 
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FIG. 2. Loci of the fixed points of Eq. (31) in the (d, a) plane, 
for b — 3.5, 5, and 10. Dashed lines: unstable branches. 



(resp. low-) a* branches, corresponding to stable steady 
growth modes, present end points for D/X = d max (resp. 
rfmin). If D/X > dmaxj the only stable growth mode is 
the a « 1, slow one; analogously, for D/X < d min , the 
only steady growth mode is the a ~ 1, fast one. Between 
these two values, depending on the initial conditions, the 
system settles in either of these two regimes. 

As discussed at the end of Section II, when tempera- 
ture T increases, b, hence <i max , decrease, while due to the 
decoupling effect A [sec Eq. (26)] decreases: A <i max thus 
decreases with T. Moreover, as shown in Ref. 9 , D ex- 
hibits a weak increase in the lOOnm range. Consequently, 
the fast growth mode exists only up to a maximum tem- 
perature T max such that: 



D 
X 



c(^max) 



(34) 



T=T n , 



Similarly, steady slow growth exists only above the tem- 
perature T m j n < T max such that: D/X = d m - ln . 

Our rugged spherulitc model thus predicts the exis- 
tence of a finite temperature range where both slow or 
fast growth can be observed. Therefore, we expect the 
following hysteretic behaviour: after a quench at very low 
temperature, a spherulite displays steady fast growth; 
as temperature is slowly increased, the front velocity re- 
mains quasi-constant before reaching the vicinity of T max 
where it slightly decreases; at T = T max , growth com- 
mutes to the slow (normal) mode 29 . Conversely, a slowly 
growing spherulite should commute to fast growth below 
T = T min < T max . The OTP data on Fig. 1 of Ref. 20 
indeed seem to present a finite (narrow) range of tem- 
peratures where both growth modes are observed. 

The value of the fast-to-slow velocity ratio Ut/u a ~ 
cosh(6) is the same as for the ideal spherulite. As seen in 



Section II, the experimental value Uf/u s ~ 10 4 for OTP, 
leads to b « 10 consistent with a reasonable estimate for 
the activation volume ao- We now examine whether our 
model is compatible with the experimental observation 
that the fast-to-slow transition occurs close to T g . For 



b = 10, we find d ma5 

at T g for OTP (V Q ~ 10" 12 m/s, 17/ju ~ 60s) this leads to 
D = Vq d max r)/[i ~ lOOnm, a quite satisfactory estimate. 



1500; with the above estimates 
-12„ 



IV. DISCUSSION 

On the basis of the above analysis, we propose that GC 
growth is a stress-induced phenomenon. Our argument 
can be summarized as follows. 

(i) Since the crystal phase is denser than the amorphous 
one, the crystallization process generates stresses, in par- 
ticular along the growth front. The growing spherulite 
can be viewed as an Eshelby inclusion, and as such, due to 
its globally spherical shape, docs not create any net free- 
volume change in the matrix. However, a shear stress 
E alone is able to accelerate interfacial kinetics, with 
respect to the normal diffusion-controlled regime, by a 
factor which we roughly estimate as ~ cosh(Sao/fcsT), 
with ao an Eyring-like activation volume, which must be 
of order a molecular volume. With interfacial stresses of 
order of the Eshelby stress for a sphere £0 = f/^Po/poj 
it immediately appears that moderate volume contrac- 
tions of order a few % are sufficient to amplify growth 
velocity by several decades (~ 10 4 in the case of OTP). 
Moreover, this relation between the level of Eshelby inter- 
facial shear stresses generated by the phase change and 
the amplification factor values is fully consistent with the 
observation by the Madison group 12 that the velocity of 
GC growth compares with that of solid-solid transforma- 
tions between polymorphs of the same material. 

(ii) Since the matrix is visco-clastic, it is able to partially 
relax the stresses generated by crystal growth. From this 
viewpoint slow growth is possible only as long as stress 
relaxation is much faster than stress generation. As dis- 
cussed in the previous section, ahead of the growth front, 
stresses decrease on a characteristic length scale £ which 
is fixed by the local geometry of the crystal. The charac- 
teristic time scale of stress variation in the matrix ahead 
of a front moving at velocity Vq is thus r ~ £/Vq. Slow 
growth, which is the rule above ~ 1.2 T g , can hence be 
sustained so long as this time scale is larger than r a , i.e.: 



V r a < £ 



(35) 



Note that since Vq and r Q scale respectively as D and 77, 
VoT a ~ rjD. For strong glassformers, the Stokcs-Einstein 
scaling 77 D ~ T holds; Vot q thus decreases with temper- 
ature: commutation from slow to fast growth cannot oc- 
cur upon cooling. It is only in fragile glassformers, where 
below T ~ 1.2 T g , rjD ~ t? 1 ^ (with < £ < 1) increases 
rapidly upon cooling, that the emergence of GC growth 
is possible. 

(iii) The last question is then to decide on the order of 
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magnitude of the length scale I. We have shown in Sec- 
tion II that modelling a spherulite as a compact homo- 
geneous sphere of radius R is too naive: since it yields 
I ~ R, GC growth can only be a transient phenomenon; 
moreover, close to T g , where GC growth is found to 
emerge, the maximum radius of a fast spherulite would 
then lie in the sub-micrometric range, contrary to obser- 
vations. We have then argued that the relevant length 
scale is the width D of the elementary substructures pro- 
duced by the intermittent branching growth process, spe- 
cific of sphcrulitic growth. We have then shown in Sec- 
tion III that the order of magnitude D ~ lOOnm is con- 
sistent with the observed value of the emergence temper- 
ature T t —T g . 

On this basis, we claim that two conditions are nec- 
essary for the emergence of steady GC growth: (i) the 
existence of the decoupling phenomenon^' 5 in the under- 
cooled liquid (i.e. a large fragility index); (ii) a fine sub- 
structure of the growing crystal (i.e. a rugged interface), 
which is the case of spherulites 9 . 

At this stage, our model is consistent with the main ex- 
perimental facts. Moreover, one of its outcomes, namely 
the prediction that the slow to GC growth "transition" 
is hysteretic, should provide an independent test of its 
validity. Experimental investigation of this question, if 
feasible, would be highly desirable. 

As a last remark, let us emphasize that our rugged 
spherulite model remains oversimplified, in particular 
as it assumes that the interfacial sub-structures share 
a common length scale D. In practice, in an actual 
spherulite, the size of sub-structures and growth veloc- 
ities of the corresponding protrusions fluctuate. Conse- 
quently, the temperature at which GC commutes to nor- 
mal growth may fluctuate along the growth front. This 
might explain the important observation 12,20 that at tem- 
peratures between T t and ~ 1.2 T g , sparse fast growing 
fibers emerge from the slow growth fronts with which 
they coexist. 



Appendix A: Appendix 

We recall here a few results about the classical Eshelby 
problem. Within an infinite elastically isotropic medium, 
we consider a domain C (the inclusion) that undergoes a 
density change as it transforms into another phase with 
the same elastic constants, under the condition of conti- 
nuity of displacement at the interface. Since the trans- 
formation is a pure compression the total resulting strain 
field can be written as 14,15 : 



12tt(1 - v) 



VV0 



(Al) 



with e = the relative density change and with <j> the 
elastic potential defined as 



1 



(A2) 



which verifies A<f> = 4 n in the inclusion and A(f> = 
outside. 

This last equation entails that V-e = in the sur- 
rounding matrix. That is, whatever the shape of the in- 
clusion, the transformation generates neither free- volume 
nor pressure field. 

Moreover, given a partitioning of the inclusion, it 
immediately results from the form of equations (Al) 
and (A2) that the strain e can be written as the sum 
of the contributions of different sub-domains. In the case 
illustrated on Fig. 3 of an inclusion formed of a sphere 
S of radius R decorated on its surface by a bulge £>, the 
stress in the matrix is thus written as a = 2p (eg + eg) 
with Is and eg the strains generated by the sphere and 
bulge respectively. 

We note that for a inclusion with rotational symmetry 
about the z axis, and defined in cylindrical coordinates 
by p < pc{z) with z € [z m ini z max ]) the zz component of 
the stress at a point f= (0, 0, z) reads: 



pel 

T i 



dz' 



(p< 2 + (z- z'f) 



3/2 



Pc(z') 



(A3) 



J 


i 






s \ 








► 
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FIG. 3. A spherical inclusion with a small hemispherical 
bulge. 

We now specialize to the case of a hemispherical bulge 
of radius D <C R, in which case we can consider the inclu- 
sion as combining a full sphere with a hemisphere. For a 
sphere of radius R centered at the origin one immediately 
recover from (A3) the classical result: 



for \z\ > R 



with cr° = the value of the tensile stress at the 

9 1— v 

surface of the sphere. For a hemisphere of radius D cen- 
tered at z = —R, we find: 



D 



z + R\ 



for z<-R-D 



with the prefactor c = 1 — ~ 0.12. At the interface, 
z = ~R—D, the tensile stress is thus ((l+D/R)~ 3 +c) a . 
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The presence of a small bulge hence only slightly mod- 
ifies the interfacial stress as compared with the value cr° 
for a perfect sphere. However, close to the interface (up 
to distances of order D at least), the stress gradient is 
dominated by the contribution of the bulge, that is, the 
length scale of the stress variations is the bulge radius D. 
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